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摘要 : 【目的 ] 旨 在 了 解 烟 蚜 Myzus persicae 不 同 翅 型 (有 翅 和 无 翅 ) 的 转录 组 特征 ,丰富 烟 蚜 转录 组 
数据 信息 ,为 探 明 烟 蚜 翅 型 分 化 的 分 子 机 制 打下 基础 。【 方 法 】 利用 高 通 量 测序 技术 Humina 
HiSeq 2000 分 别 检测 烟 蚜 车 蚜 、 无 起 成 蚜 和 有 起 成 蚜 转 录 组 表达 差异 ;通过 实时 荧光 定量 PCR 
(qPCR ) 方 法 ,分 析 翅 型 相关 基因 表达 模式 。【 结 果 】 经 过 测序 获得 质量 值 高 于 20 的 碱 基 上 比例 
(Q20) 均 不 低 于 94% 的 数据 (GenBank 登录 号 : SRR6112438, SRR6112436 和 SRR6112437) ; 所 有 
reads 组 装 成 128 065 个 unigenes ,平均 长 度 607.08 bp; 将 所 得 序列 注释 到 NR，NT,KOG，GO 和 
KEGG 等 数据 库 进行 比 对 ,共有 128 065 个 unigenes 被 注释 。 基 因 表 达 分 析 发 现 , 与 无 翅 成 蚜 和 若 
是 两 者 相 比 ,有 起 成 是 直 形成、 肽 类 激素 和 受 体 蛋 白 相关 基因 的 表达 量 存在 显著 差异 ,表明 这 些 差 
异 表达 基因 在 烟 蚜 生长 发 育 过 程 中 和 烟 蚜 如 的 形成 有 关 。qPCR 分 析 结 果 显示 ,在 有 起 成 是 与 无 起 
成 蚜 中 ,Bursicon ( bur-a 和 bur-8) , flightin (Flin), Wntless (Wnt), Distal-less ( DIL) , Decapentaplegic 
(Dpp) , Juvenile hormone epoxide hydrolase 基因 (JHEH) 和 Ecdysone receptor 基因 (kcR) 表达 量 存 在 
显著 差异 ,表明 这 些 基因 可 能 与 起 型 发 育 有 关 。[【 结论] 烟 蚜 起 形成 的 过 程 中 体内 bur-a 和 bur-B, 
Fin, Wnt, Dll, Dpp, JHEH 和 EcR 基因 表达 发 生 了 变化 。 该 研究 结果 揭示 了 烟 蚜 起 形成 相关 基因 
的 总 体 表达 特征 ,为 进一步 深入 开展 烟 蚜 起 形 成 的 研究 英 定 基础 
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Abstract: [Aim] This study aims to understand the transcriptome features of Myzus persicae of different 
wing types ( winged and apterous ), to enrich its transcriptome data, and to reveal the molecular 


mechanisms of wing dimorphism of this insect. [Methods] The differences in transcriptome expression 
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between nymphs, apterous adults, and winged adults of M. persicae were detected by high-throughput 
sequencing technology Illumina HiSeq'" 2000. The expression patterns of wing type related genes were 
analyzed by real-time quantitative PCR (qPCR). [Results] The percentage of nucleotides with the 
quality score over 20 was not less than 94% ( GenBank accession numbers; SRR6112438 , SRR6112436 
and SRR6112437). All reads were assembled into 128 065 unigenes with an average length of 607. 08 
bp. The resulting sequences were annotated to databases NR, NT, KOG, GO, KEGG, etc., and a total 
of 128 065 unigenes were annotated. Gene expression analysis showed that the genes of flight proteins, 
peptide hormones, receptor proteins, etc., which are related to wing formation during growth, were 
differentially expressed between winged adults and apterous adults or nymphs. qPCR analysis showed that 
there were significant differences in the expression levels of wing type differentiation related genes 
including Bursicon (bur-a and burB), flightin ( Fln), Wntless ( Wnt), ( DIL) , 
Decapentaplegic ( Dpp) , juvenile hormone epoxide hydrolase gene ( JHEH) and ecdysone receptor gene 


Distal-less 


(EcR) between winged adults and apterous adults. [ Conclusion] The expression levels of bur-a, bur-B, 
Fin, Wnt, DIL, Dpp, JHEH and EcR change during wing formation in M. persicae. Our results show the 
general expression characteristics of wing formation genes and provide a basis for further research of wing 
formation in M. persicae. 


Key words: Myzus persicae; transcriptome; wing development; differentially expressed gene; high- 


throughput sequencing 





烟 蚜 Myzus persicae 又 名 桃 是 ,是 一 种 危害 性 很 
大 的 农业 害虫 , 它 寄主 广泛 ,除了 危害 烟草 ,同时 也 
危害 多 种 农作物 。 烟 蚜 不 仅 利 用 刺 吸 式 口 絮 吸取 烟 
株 汁液 对 烟 株 造成 直接 危害 外 ,还 可 传播 多 种 烟草 
病毒 病 ,间接 对 烟草 造成 更 严重 的 危害 ( 秦 焕 菊 和 
EHEZ, 1996)。 目 前 ,在 田间 对 烟 蚜 的 综合 治理 以 
化 学 防治 为 主 , 但 杀 虫 剂 不 合理 使 用 和 施 药 器 械 落 
后 等 原因 使 得 烟 蚜 对 杀 虫 剂 长 期 处 于 高 选择 和 高 残 
留 环 境 中 ,加 上 烟 蚜 具有 较 强 的 环境 适应 能 力 和 较 短 
的 生活 史 , 抗 ( 耐 ) 药 性 问题 亦 日 趋 严重 (Devonshire et 
al., 1998) , 烟 蚜 的 防治 问题 也 日 趋 严 峻 。 因 而 , 寻 
找 新 型 的 防治 技术 防治 烟 蚜 势 在 必 行 。 

烟 蚜 是 迁 飞 性 昆虫 , 当 周 围 环境 不 利于 烟 蚜 的 
生存 时 ,就 会 形成 有 示 型 迁 往 适 宜 的 生存 环境 ( 刘 
向 东 等 , 2004) 。 无 翅 烟 师 和 有 翅 烟 是 作为 两 个 不 
同形 态 的 个 体 ,其 生长 发 育 时 期 体内 的 生理 状况 也 
不 一 样 ( 王 凯 等 , 2008) 。 在 烟 蚜 在 不 同 发 育 时 期 和 
有 却 型 烟 蚜 的 发 育 过 程 中 ,其 体内 相关 基因 表达 也 
会 相应 发 生 了 变化 ,但 对 其 差异 基因 表达 及 相关 调 
控 机 理 尚 未 进行 深入 人 研究。 

近年 来 , 随 着 高 通 量 测序 技术 和 转录 组 学 的 发 
展 与 不 断 完 善 〈 张 棋 麟 和 彭 明 龙 , 2013; 胡 俊 杰 等 ， 
2017) ,蚜虫 转录 组 研究 取得 了 一 定 的 研究 成 果 
( Costechareyre et al., 2013; Anathakrishnan et al., 


2014)。 但 关于 烟 蚜 转录 组 研究 报道 较 少 ,对 烟 蚜 









































不 同 翅 型 的 转录 组 表达 信息 的 了 解 甚 少 ,尤其 是 却 
的 形成 及 发 育 阶 段 。 本 研究 拟 通过 转录 组 测序 , 获 
得 烟 蚜 在 不 同 虫 态 及 翅 型 的 转录 组 基本 信息 ,并 初 
步 分 析 烟 蚜 若 蚜 .无 翅 成 蚜 和 有 却 成 是 的 转录 组 信 
息 ,了 解 烟 蚜 在 生长 发 育 过 程 中 特定 基因 的 表达 情 
DL ,为 深入 研究 烟 师 翅 型 分 化 机 制 和 相关 基因 克隆 
及 功能 分 析 葛 定 一 定 的 基础 。 


























1 材料 与 方法 
1.1 供 试 虫 源 及 饲养 方法 
供 试 虫 源 : 烟 蚜 于 2008 年 采 自 贵州 省 长 顺 县 马 
路 乡 烟叶 生产 基地 (26°03'N, 106?45'E; 海拔 1 280 
m) ,在 不 接触 任何 药剂 的 情况 下 利用 烟 株 连续 饲 
养 了 9 年 。 饲 养 条 件 : 温度 25 + 1% ,相对 湿度 
60% ~70% , 光 周 期 14L: 10D。 
1.2 主要 试剂 和 仪器 

Total RNA Extractor ( Trizol) ( Sangon Biotech, 
上 海 ), Qubit2. 0 RNA 483 it] x (Life, 上 海 )， 
VAHTS" mRNA-seq V2 Library Prep Kit for 
Illumina 9 和 VAHTS"" DNA Clean Beads ( Vazyme 
Biotech ,南京 ) 。 

Qubit2. 0 2€ 363] ( Invitrogen ,美国 ) ,微型 洲 涡 
混合 仪 (上 上海 沪 西 分 析 仪器 厂 ) ,台式 高 速 低温 离心 
机 (Thermo ,美国 ) ,PCR 仪 (BIO-RAD ,美国 ) ,电泳 
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仪 (北京 市 六 一 仪器 厂 ) 。 
1.3 样品 的 采集 

AIA (nymph) 的 采集 :将 同一 批 初 产 烟 蚜 奎 蚜 
(1-4 龄 各 25 头 ) 置 于 一 干净 的 离心 管内 , 放 入 
-80C 冰箱 中 保存 , 待 提取 RNA。 无 翅 成 蚜 
(apterous adult) 的 采集 :将 同一 批 4 ti ep ELT S 
一 株 无 虫 的 健康 烟 株 上 ,将 烟 株 置 于 适宜 饲养 条 件 
(温度 : 25 +1%C 、 相 对 湿度 :60% ~ 7096 、 光 周期 
14L: 10D ) 下 待 若 蚜 虹 皮 形成 成 蚜 , 挑 取 25 头 1 日 
龄 无 起 成 蚜 置 于 干净 的 离心 管内 , 放 入 - 80Y 冰箱 
中 保存 , 待 提取 RNA, ARIY (winged adult) 的 采 
集 :将 同一 批 4 龄 车 蚜 置 于 男 一 株 无 虫 的 健康 烟 株 
上 ,将 烟 株 置 于 适宜 饲养 条 件 ( 温度 25 x 1C 、 相 对 
湿度 60% ~70% 、 光 周期 14L: 10D) FP frt Bp Ez 
形成 成 蚜 , 挑 取 25 头 工 日 龄 有 翅 成 蚜 置 于 干净 的 离 
心 管内 , 放 入 -80 冰箱 中 保存 , 待 提取 RNA, 
1.4 样品 RNA 的 提取 ( Trizol 法 ) 

将 1.3 节 3 组 样品 按照 RNA 提取 的 标准 方法 
分 别提 取 烟 是 若是 .无 翅 成 师 和 有 翅 成 师 样 品 的 总 
RNA ,并 进行 RNA 质量 检测 。 合 格 后 ,将 每 组 置 于 
- 80C 冰箱 中 保存 待 用 。 
1.5 cDNA 文库 的 建立 和 转录 组 测序 

利用 Qubit2.0 RNA 检测 试剂 盒 对 总 RNA 精确 
定量 ,以 确定 文库 构建 所 加 入 总 RNA 的 量 ,利用 
Qubi2. 0 DNA 检测 试剂 盒 对 回收 的 DNA 精确 定 
量 , 以 便 按照 1:1 的 比例 等 量 混合 后 在 Humina 
HiSeq “2000 测序 平台 上 进行 测序 。 并 对 测序 获得 
原始 数据 经 过 去 除 低 质量 片段 ,切除 reads 中 含 N 
序列 后 ,获取 纯净 序列 。 具 体 参 照 Wang 等 (2012 ) 
研究 方法 。 
1.6 不 同 虫 态 烟 蚜 差异 基因 表达 分 析 

将 试验 获取 的 Clean reads 转录 组 数据 进行 混 
合 拼接 ,然后 利用 FPKM 法 计算 每 个 unigene 的 表达 
量 (Mortazavi et al., 2008; 吴 小 波 等 , 2014), 

在 此 基础 上 ,参照 Audic 和 Claverie( 1997 ) 基于 
测序 的 差异 基因 检测 方法 ,对 若是 无 翅 成 蚜 和 有 却 
成 蚜 差 异 表达 基因 进行 躲 选 。 最 后 采用 FDR 的 方 
法 对 差异 检验 的 己 值 作 多 重 假设 检验 校正 得 到 0 
值 ,为 了 得 到 显著 差异 表达 的 基因 ,我 们 将 筛选 条 件 
设 为 :0 值 <0.001 且 差 异 倍数 1Fold changel >2。 分 
析 这 些 差异 基因 在 GO (Gene Ontology) 中 的 分 布 情 
况 并 阐明 样本 差异 在 基因 功能 上 的 体现 ,在 KEGG 
(Kyoto Encyclopedia of Genes and Genomes ) 数据 库 分 
析 差 异 表达 基因 参与 的 最 主要 生化 代谢 途径 和 信号 









































转 导 途 径 。 
1.7 差异 表达 基因 的 荧光 定量 PCR 验证 

为 了 验证 试验 结果 的 可 靠 性 ,根据 相关 文献 推测 
有 翅 型 和 无 翅 型 烟 蚜 中 差异 表达 基因 与 怒 发 育 的 关 
Z ,挑选 与 翅 生 长 发 育 关 系 较 密切 的 基因 进行 荧光 定 
量 PCR 检测 。 按 照 1.3 节 样 品 采集 方法 分 别 重 新 采 
集 样 品 ( 每 个 样品 设 3 次 生物 学 重复 和 3 KERE 
复 ) ,按照 1.4 节 样 品 制备 的 方法 (Trizol 法 ) 制备 
qPCR 验证 试验 样品 。 从 转录 组 数据 中 挑选 8 个 有 关 
飞行 和 怒 型 发 育 的 基因 ,分 别 为 Bursicon ( bur-a 和 
bur-B), flightin ( Fin), Wniless ( Wnt), Distal-less 
(Dll), Decapentaplegic ( Dpp ), Juvenile hormone 
epoxide hydrolase 基因 (JHEH) 和 Ecdysone receptor 基 
(EcR) ,以 B-actin ( GenBank 登录 号 : MF083568) 为 
内 参 基因 进行 荧光 定量 PCR, PCR 反应 体系 :iQ™ 
SYBR 9 Green Supermix 10 pL, 正 反 向 引物 (10 
umol/L) 各 1 uL, cDNA 模板 1 pL, ddH,0 7 AL。 
PCR 反应 程序 :95%C 预 变 性 10 min; 95% 变性 30 s, 
55% 退火 并 延伸 30 s, 40 个 循环 ;接着 65 ~95%C 的 溶 
解 曲线 。 检 测 分 析 它 们 在 若 蚜 、 无 翅 成 蚜 和 有 起 成 蚜 
中 的 相对 表达 量 ,最 后 利用 2 "法 (Livak and 
Schmittgen，2001 ) 计算 相对 表达 量 ,以 探讨 这 些 基 
与 茵 发 育 的 关系 。 引 物 设计 如 表 1。 






































表 1 实时 荧光 定量 PCR 引物 
Table 1 qPCR primers 





引物 引物 序列 (5 -3') 
Primers Primer sequences 
yg-Fln-F AACCAGTTGAAGCTCCTAAACC 

yg -Fln-R AGTCGTCGTAGTAATTGCTCTG 
yg-EcR-F ACCTCGTCGTCATCTTCCTC 
yg-EcR-R GCTGCTGATGGTTGGTGTAG 
yg-Wnt-F CTGCTGTCGTCGGTGGTAAT 
yg-Wnt-R TAGGTGACTGTTCCAGGTAGAC 
yg-Dpp-F CGGCTACGACGCTTACTACT 
yg-Dpp-R GCATAGATATGGACGCCAACTG 
yg-DII-F TAAGCGGTGTGCTGTTGATG 
yg-DII-R CGGCTGGTATCTCCACTTTG 
yg-Eh-F AACATCACGACCGAAGTCATC 
yg-Eh-R CGAGTTCAGGTACGCTTGTC 
yg-burs-oF GTCATCCATGTTCTACAATAC 
yg-burs-oR CTTCTCGTTCACCACTCTC 
yg-burs-BF TGGTGACATAGCAGTGAAC 
yg-burs-BR GGATCGTAGCAATGTGTAAG 
Actin-F TGGTATCGTCTTGGATTCTG 
Actin-R TTAGGTAGTCGGTGAGATCA 














2 Hj 
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1.8 数据 统计 分 析 
数据 采用 Microsoft Excel 2007 进行 数据 处 理 ,再 
用 SPSS16.0 软件 进行 单 因素 方差 分 析 和 多 重 比较 。 





2 结果 


2.1 烟 蚜 转录 组 Illumina 测序 和 序列 组 装 
对 烟 蚜 若是 、 无 翅 成 是 和 有 起 成 蚜 样 品 进 行 
Illumina HiSeq 2000 测序 ,将 测序 获得 的 原始 数据 


(GenBank 登录 号 : SRR6112438, SRR6112436 和 
SRR6112437 ) 文件 经 CASAVA 碱 基 识 别 (base 
calling) 分 析 转 化 为 原始 测序 序列 (sequenced 
reads) ,我 们 称 之 为 raw data 或 raw reads, LJ FASTQ 
(简称 为 fq) 文件 格式 存储 , 其 中 包含 测序 序列 
(reads ) 的 序列 信息 以 及 其 对 应 的 测序 质量 信息 ;并 
使 用 Trimmomatic 进行 数据 处 理 , 过 滤 去 杂 ( 包 括 含 
有 带 接头 .重复 的 和 质量 很 低 的 序列 ) ,得 到 纯净 的 
序列 ( 表 2) 。 














表 2， 烟 蚜 转录 组 数据 分 析 
Table2 Analysis of the transcriptome data of Myzus persicae 
"S 所 有 reads 数 所 有 碱 基数 PAFI PAZERA bp) N 碱 基 比 例 (% ) Xm 
样本 GC 含量 (%) 
(bp) (bp) (bp) Total clean Q20 Proportion of 
Sample GC content 
Total reads Total base count Total clean reads nucleotides N bases 
有 却 成 蚜 Winged adult — 48 594 342 7 289 151 300 48 217 182 6 905 539 640 94. 003 0.002 47.485 
若 蚜 Nymph 64 720 056 9 708 008 400 64 331 519 9 229 450 304 94.428 0. 002 42.525 
JERES Apterous adult — 65 466 494 9 819 974 100 65 099 307 9 425 647 502 94.891 0. 005 37.177 


Q20: 碱 基质 量 在 20 以 上 的 比例 Proportion of bases with the quality score over 20. 


从 表 2 中 可 以 看 出 ,经 过 测序 获得 总 序列 数 为 
48 594 342 ~65 466 494 bp 之 间 , 其 中 干净 序列 数 在 
48 217 182 ~65 099 307 bp 之 间 , 干 净 序 列 长 度 介 于 
6 905 539 640 ~9 425 647 502 bp 之 间 。 碱 基质 量 
1E 20 UE (Q20) 的 碱 基 比例 均 高 于 94% ,GC 含量 
比例 为 37.177% ~47.485% 之 间 ( 表 2)。 对 测序 数 
据 进行 组 装 后 分 析 表 明 ( 表 3 ) ,本 研究 获取 的 样品 





128 065 个 unigenes ( 表 4)。 其 中 注释 到 KOG, 
NR, NT, GO 和 KEGG 等 数据 库 的 结果 分 别 是 
27 755 (21. 6796), 48 403 (37. 896), 48 359 
(37.7690), 40 556 (31. 6796) 和 19 995 
(15.6196) , f Hj NCBI Blast 将 unigene 基因 序列 
E NR 数据 库 进 行 同 源 性 比 对 , 取 耻 vlaue « 1e 7I 
相似 度 >30% 的 比 对 结果 ,得 到 unigene 的 功能 注释 























所 有 序列 能 组 装 成 145 200 条 Transcript 转录 本 序 
列 ,平均 长 度 742. 45 bp, Transcript 的 N50 为 1 477 
bp。 将 所 得 样品 的 clean reads 数据 de novo 利用 
Trinity 进行 组 装 总 共 获 得 128 065 条 unigenes 转录 
本 序列 ,unigenes 的 N50 为 855 bp, 3£35] KK BE 607. 08 
bp( 表 3 ) 。 对 最 后 组 装 出 来 的 unigene 进行 长 度 分 布 
统计 ,长 度 超过 了 1 000 bp 共有 20 287 条 (图 1)。 
2.2 Unigene 功能 注释 

将 unigene 分 别 注释 到 KOG, NR, NT, GO 和 
KEGG 等 数据 库 中 ,并 分 别 对 注释 到 每 个 库 以 及 注 
释 上 的 unigene 数目 进行 统计 。 本 实验 总 共 获 得 








言 息 和 同 源 物种 信息 (图 2)。 结 果 显示 ,分 别 有 

24. 1396 , 23. 6596 和 9. 06% 的 基因 与 小 家 鼠 Mus 
musculus , Wi =. Sf Acyrthosiphon pisum FI 2e XX FÆ f 
Diuraphis noxia 比 对 上 ， 另外 还 有 23. 93% 的 未 知 基 
因 ( 图 2)。 

将 这 些 unigenes 在 KOG 数据 库 中 进行 分 类 , 共 
分 26 类 (图 3) ,其 中 最 大 群体 是 R 类 (一 般 功能 预 
测 ) ,有 3 641 个 基因 ,其 次 是 T, 0, ] 入 类 (分 别 
为 “信号 转 导 机 制 " “蛋白质 翻转 和 分 子 伴侣 类 ”， 
“翻译 .核糖 体 结构 与 生物 起 源 " 和 ”转录 ”) ,分别 对 
应 3 513, 2 644, 1 876 和 1 855 个 基因 (图 3)。 


























RI 烟 蚜 转 录 组 测序 序列 组 装 质量 
Table 3 Assembly quality of Myzus persicae transcriptome 





长 度 范 围 





最 大 长 度 (bp) 最 小 长 度 (bp) 








总 数 N50 长 度 (bp) 2 总 长 度 (bp) ”平均 长 度 (bp) 
Length range Maximam Minimum 
Total number N50 length Total length Mean length 
z500bp zl 000 bp length length 
转录 本 Transcript 145 200 46 464 23 317 1477 118 446 201 107 804 096 742.45 
Unigene 128 065 33 844 14 164 855 118 446 201 77 745 381 607.08 
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图 1 Transcript 和 unigene 长 度 分 布 
Fig. 1 Length distribution of transcripts and unigenes 
表 4 烟 蚜 基因 注释 结果 统计 
Table 4 Annotation of unigenes of Myzus persicae 
数据 库 注释 的 unigenes 数目 注释 的 unigenes 百分比 
Database Number of annotated unigenes Percentage of annotated unigenes 
CDD 32 823 25.63 
KOG 27 755 21.67 
NR 48 403 37.8 
NT 48 359 37.76 
PFAM 28 714 22.42 
Swiss-Prot 39 320 30.7 
TrEMBL 48 953 38.23 
GO 40 556 31.67 
KEGG 19 995 15.61 
在 至 少 一 个 数据 库 中 注释 成 功 
64 766 50.57 
Annotated in at least one database 
在 所 有 数据 库 中 都 注释 成 功 
12 034 9.4 
Annotated in all database 
总 基因 数目 Total unigenes 128 065 100 














Ar T AEA 无 翅 成 师 和 有 翅 成 师 中 差异 表达 基 
,具体 差异 数量 如 图 5 所 示 , 从 图 中 可 以 看 出 ,有 
翅 成 蚜 中 上 调 表达 的 基因 数目 均 比 奋 虹 和 无 翅 成 是 
中 多 ;无 翅 成 蚜 与 若是 相 比 ,无 翅 成 虹 中 上 调 表达 的 
基因 数目 要 比 知 蚜 中 少 。 


进一步 将 注释 到 GO 数据 库 中 的 394 077 个 
unigenes 分 成 生物 学 过 程 细胞 组 分 、 分 子 功 能 三 大 
类 别 , 其 中 生物 学 过 程 分 成 23 个 类 别 ,细胞 组 分 分 
成 19 个 类 别 ,分 子 功能 分 成 20 个 类 别 ( 图 4)。 

2.3 烟 蚜 不 同 翅 型 中 的 差异 表达 基因 
为 了 检测 烟 师 不 同 虫 态 基因 表达 差异 ,我 们 分 

















张 ” 超 等 : MFE ANKA SUMI JOH DRUUEEG S62 BORDURE RED B] 36 35 25 54) T 161 





INX Mus musculus 

Bii 28 Acyrthosiphon pisum 

麦 双 尾 蚜 Diuraphis noxia 

人 类 Homo sapiens 

WIRE Diachasma alloeum 

OSF Pantholops hodgsonii 
WX Rattus norvegicus 

Ea] E PRESE Fopius arisanus 
NAE Microplitis demolitor 
IKE Cricetulus griseus 

TEW Lasius niger 

RS Columba livia 

ZZE Cephus cinctus 
TiSIDEZXIB. Trichuris trichiura 
HIR Plutella xylostella 


RIE Macaca fascicularis 
EEDS Peromyscus maniculatus bairdii 






24.13% 











9.06% 























23.93% 4 & BK Mesocricetus auratus 
2 Nematostella vectensis 


其 他 Others 


2.7896 


图 2 NR 数据 库 同 源 物 种 分 类 图 


Fig. 2 Species distribution of homology search of unigenes against the NR database 


ARNA 加 工 和 修饰 RNA processing and modification 

B 染色 质 结构 和 动力 学 Chromatin structure and dynamics 

C 能 量 生 产 čt Energy production and conversion 

fuu, REA Cell cycle control, cell division, chromosome partitioning 
谢 Amino acid transport and metabolism 

谢 Nucleotide transport and metabolism 

与 代谢 Carbohydrate transport and metabolism 

H 4s 谢 Coenzyme transport and metabolism 





















1 脂 类 转运 和 代谢 Lipid transport and metabolism 
J 翻译 、 核 糖 体 结构 与 牛 物 起 源 Translation, ribosomal structure and biogenesis 
K 转录 Transcription 
工 复 制 ， 重 组 和 修复 Replication, recombination and repair 
M 细胞 壁 / 膜 / 包 膜 牛 物 起 源 Cell wall/membrane/envelope biogenesis 
细胞 Cell motility 
口 翻 译 后 修饰 、 贷 白质 翻转 和 分 子 伴 但 类 Posttranslational modification, protein turnover, chaperones 
P 了 无 机 离 了 转运 与 代谢 Inorganic ion transport and metabolism 
Q 次 生 代谢 产物 的 生物 合成 、 运 输 和 分 解 代谢 Secondary metabolites biosynthesis, transport and catabolism 
R 一 般 功 能 预测 General function prediction only 









Hi ATWATDAE TRIS Intracellular trafficking, secretion, and vesicular transport 
V 防御 机 制 Defense mechanisms 

W 胞 外 结构 Extracellular structures 

X AAA Unnamed protein 

立 核 结构 Nuclear structure 

Z 细胞 骨架 Cytoskeleton 
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2776 


基因 数量 Number of genes 


基因 所 占 百 分 比 Percentage of genes 


3 unigenes 的 KOG 分 类 图 
Fig. 3 Clusters of orthologous groups ( KOG) classification of the unigenes 
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unigenes 的 GO 分 布 图 


图 4 


Fig. 4 Gene ontology classification of the unigenes 
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100007 m -3H Up 上 调 表达 ,3 197 个 基因 下 调 表达 ,包括 66 个 有 关 翅 
3 SA 型 生长 发 育 的 基因 (其 中 40 个 上 调 ,26 个 下 调 ) ( 表 
S nd 6) 和 35 个 有 关 表皮 生长 的 基因 (其 中 26 个 上 调 ,9 
i FN 个 下 调 ) 。 若 蚜 与 无 翅 成 蚜 相 比 , 它 们 之 间 有 10 829 
z 个 基因 存在 差异 表达 ,其 中 在 若 蚜 中 上 调 表达 的 基 
Ri 因 有 6 422 个 ,在 若 蚜 中 下 调 表达 的 基因 有 4 407 个 
H 4000 - GET. 
E 2.4 ”差异 表达 基因 的 荧光 定量 PCR 验证 
AUN 从 差异 表达 基因 中 选择 了 8 个 关于 飞行 和 起 发 
Y 育 的 基因 在 若 蚜 Do JUBERUGS ARAE m rh EG 
Wos N Wos A Nus A qPCR 验证 。 实 验 结 果 表 明 , Bursicon 基因 (bxr-a 和 
图 5 烟 是 若是 ,有 址 成 蚜 和 无 翅 成 蚜 间 的 bur-B) , Flightin 基因 (Fin)，Wntless 基因 ( Wnt)， 
差异 表达 基因 数目 Distal-less 基因 ( DII) , Decapentaplegic 基因 ( Dpp) , 


Fig. 5 Numbers of differentially expressed genes ( DEGs) 
in nymphs, winged adults and apterous adults 
of Myzus Persicae 
W vs N: ARIS vs #1 Winged adult vs nymph; W vs A; A% 
成 蚜 vs 无 翅 成 蚜 Winged adult vs apterous adult; N vs A; 若 蚜 vs 
无 翅 成 蚜 Nymph vs apterous adult. 


KRA ARISE PEEL x PRI T AI I] es TE 
7 127 个 差异 基因 ,其 中 6 230 个 基因 在 有 翅 成 是 中 
上 调 ,897 个 基因 下 调 ,包括 52 个 有 关 翅 型 生长 发 
育 的 基因 (其 中 50 个 上 调 ,2 个 下 调 )( 表 5) 和 23 
个 有 关 表皮 生长 的 基因 (其 中 22 个 上 调 ,1 个 下 
调 )。 有 起 成 蚜 与 无 让 成 蚜 相 比较 ,有 12 135 个 基 
因 存 在 差异 表达 ， 其 中 8 938 个 基因 在 有 翅 成 是 中 


表 5 


Juvenile Hormone Epoxide hydrolase 基因 (JHEH) 和 
Ecdysone receptor 基因 ( EcR ) 基 在 有 起 成 蚜 体 内 
表达 量 显著 高 于 若 蚜 和 无 翅 成 蚜 中 ,其 中 Fin 基因 
表达 差异 最 为 显 善 ,有 翅 成 师 中 表达 量 是 若是 中 的 
3.34 倍 , 是 无 翅 成 蚜 中 的 95. 96 倍 ( 图 6) ,与 转录 组 
结果 一 致 。 


















































3 讨论 











本 研究 通过 高 通 量 测序 技术 对 烟 蚜 若 蚜 .无 翅 
成 师 和 有 翅 成 蚜 的 转录 组 差异 进行 了 分 析 , 获得 烟 
蚜 的 转录 组 信息 。 测 序 结果 发 现 ,所 有 试验 样本 高 
质量 的 数据 ( Q20) 均 高 于 94% ,说 明 测 序 质量 较 



































烟 蚜 转录 组 中 有 翅 成 蚜 与 若 蚜 之 间 表 达 差 异 最 显著 的 有 关 翅 型 发 育 的 12 个 基因 


Table 5 Top 12 differentially expressed genes ( DEGs) in Myzus persicae transcriptome related to 


wing development between winged adults and nymphs 


FPKM 值 FPKM value 


































































































基因 产物 基因 号 有 翅 成 蚜 vs 若 蚜 P 值 
Cene produi Gene ID AARI m Winged adieus nymph P-value 
Winged adult Nymph 

KITÆ A Flightin c120195_g1 2 254.57 30. 93 E38 U 0.001 
HÆ H C Troponin C c122884_g1 509.38 82.59 EX U 0.001 
VELLE Bursicon 25223, gl 30.61 14. 05 EJ U 0.001 
MYM 型 锌 指 蛋白 Zinc finger MYM-type c86228, gl 706.26 17.07 EJ U 0.001 
胰岛 素 受 体 Insulin receptor c26217, gl 5.04 23.71 调 Down 0.001 
章鱼 胶 受 体 B Octopamine receptor-beta c14012, gl 20.14 8.82 E U 0.001 
G 蛋白 偶 联 受 体 1 G protein-coupled receptor 1 —— 25696. gl 7.98 29.94 调 Down 0.001 
翅 极 性 蛋白 stan Wing polarity protein stan c38442, gl 4.21 20. 83 调 Down 0.001 
G 蛋白 偶 联 受 体 2 G protein-coupled receptor 2.— 70240, gl 13.77 37. 46 调 Down 0.001 
肌 球 蛋白 重 链 Myosin heavy chain c20098, gl 4.49 25.52 调 Down 0.001 
EEH Flightless c36649. gl 4.44 34.23 调 Down 0.001 
肌 球 蛋白 5Myosin 5 c37942 2 2.17 17.68 调 Down 0.001 


- : 无 差异 No difference. 3€ 6 和 7 |E] The same for Tables 6 and 7. 
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表 6 转录 组 中 有 起 成 蚜 与 无 起 成 蚜 之 间 表 达 差 异 最 显著 的 有 关 翅 型 发 育 的 12 个 基因 
Table 6 Top 12 differentially expressed genes (DEGs) in Myzus persicae transcriptome related to 


wing development between winged and apterous adults 


FPKM 值 FPKM value 




































































































































































基因 产物 基因 号 有 过 成 蚜 vs 无 区 成 蚜 P 值 
Gene product Gene ID 有 起 成 是 TURI Winged adult vs apterous adult P-value 
Winged adult Apterous adult 

KÍTÆ A Flightin c120195 gl 2 254.57 0.44 上 调 0.001 
肌 钙 蛋白 C Troponin C c122884 gl 509.38 0 调 0.001 
HALIA Bursicon c25223 gl 30.61 11.08 调 0.001 
MYM 7i £182E [1 Zinc finger MYM-type c86228 gl 706.26 4.12 调 0.001 
胰岛 素 受 体 Insulin receptor c26217. gl = 三 = 

章鱼 胺 受 体 B Octopamine receptor-beta c14012_g1 20.14 4.86 调 0.001 
G 蛋白 偶 联 受 体 1 G protein-coupled receptor 1  c25696_gl 7.98 17.02 调 Down 0.001 
翅 极 性 蛋白 stan Wing polarity protein stan c38442. gl 4.21 11.78 调 Down 0.001 
G 蛋白 偶 联 受 体 2 G protein-coupled receptor 2 —— 70240. gl = = = 

肌 球 蛋白 重 链 Myosin heavy chain c20098_gl - - 

EWEA Flightless c36649_g1 17.53 调 Down 0.001 
肌 球 蛋白 5 Myosin 5 c37942 g2 9.22 调 Down 0.001 

R7 转录 组 中 若 蚜 与 无 翅 成 蚜 之 间 表 达 差 异 最 显著 的 有 关 起 型 发 育 的 12 个 基因 
Table7 Top 12 differentially expressed genes ( DEGs) in Myzus persicae transcriptome related to 
wing development between nymphs and apterous adults 
FPKM 值 FPKM value 
基因 产物 基因 号 若 蚜 vs 无 翅 成 是 P fH 
Gene product Gene ID 无 地 成 是 Nymph w apterous adult P-value 
Nymph Apterous adult 

KÍTÆ A Flightin c120195_g1 30.93 0.44 E38 U 0.001 
肌 钙 蛋白 C Troponin C c122884 gl 82.59 0 E38 U 0.001 
IALIA Buesicon c25223 gl - - - 
MYM 型 锐 指 蛋白 Zinc finger MYM-type c86228 gl 17.07 4.12 E38 U 0.001 
胰岛 素 受 体 Insulin receptor c26217 gl 23.77 8.21 上 调 UU 0.001 
音 鱼 胺 受 体 B Octopamine receptor-beta c14012, gl - = 

G 蛋白 偶 联 受 体 1 G protein-coupled receptor 1  c25696_gl 580. 83 65.17 E U 0.001 
翅 极 性 蛋白 stan Wing polarity protein stan c38442 gl - - 

G 蛋白 偶 联 受 体 2 G protein-coupled receptor 2 — 70240. gl 37.46 17.48 E U 0.001 
肌 球 蛋白 重 链 Myosin heavy chain c20098. gl 25.52 0.53 EŠ] U 0.001 
无 翅 蛋 扬 Flightless c36649. gl E - 
肌 球 蛋白 5 Myosin 5 c37942_g2 = - 








好 。 所 有 转录 组 测序 获得 的 reads 组 装 后 , 获取 
128 065 个 unigenes ,平均 长 度 607. 08 bp , 而且 大 部 
分 unigenes 能 比 对 到 亲缘 关系 较 近 的 物种 ,这 也 说 
明 测 序 组 装 质量 较 好 。 对 有 翅 成 蚜 与 若 蚜 及 无 翅 成 
蚜 差异 表达 基因 进行 分 析 , 发 现 肽 类 激素 和 有 关 却 
的 生长 发 育 的 基因 表达 均 存 在 差异 , 且 与 翅 形 成 的 
相关 基因 显著 富 集 ,说 明 这 些 基因 与 烟 蚜 翅 的 发 育 
和 飞行 相关 。 

Flightin 作为 在 间接 飞行 肌 (indirect flight 
muscle, IFM) 中 特异 表达 的 一 种 蛋白 ,1993 年 首次 























在 果 蝇 体内 被 鉴定 出 来 ( 徐 升 胜 等 , 2009) 。 研 究 表 
BH .flightin 基因 对 于 粗 肌 丝 的 组 装 和 性 能 有 重要 意 
义 , 其 对 IFM 发 育 和 功能 中 必 不 可 少 (Vigoreaux and 
Perry, 1994; Kronert et al., 1995; Henkin et al., 
2004 ; Xue et al., 2010) ,可 见 flightin 对 昆虫 的 飞行 
有 具有 重要 作用 。 本 研究 结果 显示 ,flightin 基因 表达 
量 在 有 翅 成 蚜 中 均 高 于 无 翅 成 蚜 和 若 蚜 ,此 研究 结 
果 与 前 人 的 研究 结果 一 致 ( 张 小 琴 等 , 2013; Xue et 
al., 2013; 何 海 芳 等 ，2015 ) 。 由 此 可 见 flightin 基 
因 在 昆虫 翅 的 形成 方面 有 着 重要 的 作用 。 
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图 6 8 个 基因 在 烟 蚜 若 是 有 示 成 是 和 无 翅 成 蚜 中 的 表达 量 
Fig. 6 Expression levels of 8 genes in nymphs, winged adults and apterous adults of Myzus persicae 
Ny: 若 蚜 Nymph; Ap: 无 翅 成 蚜 Apterous adult; Wg: 有 翅 成 蚜 Winged adult. 图 中 数据 为 平均 值 上 标准 差 , 柱 上 不 同 字母 表示 差异 显著 ( 单 因素 
方差 分 机 ，LSD Duncan, P «0.05) Data in the figure are mean € SD. Bars with different superscripts differ significantly ( one-way ANOVA, LSD, 
Duncan, P «0.05). 

















E "38 AE Bi E SE HH 3 个 途径 调控 的 , 即 Kombert, 1994) , 背 腹 轴 (D-V) ,代表 基因 有 Apterous 
前 后 轴 (A-P), 人 代表 基 因 有 Engrailed (En), (Ap), Wntless ( Wnt) , Distal-less ( DIL) 等 (Cohen et 
Hedgehog (hh) , Decapentaplegic ( Dpp ) 5& ( Tabata et al., 1992) , 近 远 轴 (P-D) ,如 Aristaless 等 基因 (于 居 
al., 1992; Basler and Struhl, 1994; Tabata and JE, 2014) 。 本 实验 从 上 述 基 因 中 随机 挑选 3 个 基 
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因 进 行 验证 ,结果 显示 Wa, Dll 和 Dpp 在 有 翅 成 蚜 
中 表达 量具 高 于 无 翅 成 蚜 和 若 蚜 ,并 且 表 达 量 差异 
达到 显著 水 平 。 该 实验 结果 与 Brisson 等 (2010 ) 及 
Abouheif 和 Wray (2002) 的 研究 一 致 ,可 见 此 类 基因 
对 昆虫 码 型 的 发 育 有 着 密切 的 联系 。 

保 幼 激素 (juvenile hormone, JH) 对 昆虫 变态 及 
生殖 等 功能 具有 重要 的 调控 作用 (Riddiford, 1994; 
Wyatt and Davey，1996 ) 。 昆 虫 在 不 同 发 育 阶 段 保 
幼 激素 表达 量 的 变化 对 昆虫 变态 和 生殖 具有 必然 的 
关系 ,例如 保 幼 激 素 在 昆虫 成 虫 期 保 幼 激素 发 挥 独 
立 和 直接 的 作用 ,促进 脂肪 体 的 合成 和 生殖 发 育 
(Riddiford, 2013) o JH 调控 昆虫 的 变态 , 保 幼 激素 
环 氧 水 解 酶 (juvenile hormone epoxide hydrolase, 
JHEH) 能 降解 保 幼 激素 ,本 实验 中 有 翅 成 是 JHEH 
的 表达 量 高 于 无 翅 成 蚜 ,, 克 是 否 对 烟 蚜 翅 发 育 具有 
调控 作用 还 需 我 们 进一步 实验 验证 。 

轰 化 激素 (bursicon，pburs ) 是 调节 昆虫 表皮 硬化 
及 展翅 的 激素 类 和 蛋白 , 它 在 昆虫 表皮 硬化 和 要 伸展 
过 程 中 起 着 关键 作用 。 初 步 研 究 发 现 有 运 成 蚜 炮 化 
激素 基因 的 表达 量 均 高 于 若 蚜 和 无 运 成 蚜 , EEG 
激素 基因 受到 干扰 时 , 烟 蚜 的 这 与 表皮 不 能 正常 发 
育 ( 张 超 等 , 2017) ;与 此 同时 我 们 还 发 现 烟 蚜 蛇 皮 
激素 受 体 基因 的 表达 量 受 到 影响 。Wu 等 (2012 ) 向 
Jk Ks 4 龄 若 蚜 注射 激素 受 体 (ecdysone receptor, 
EcR) 基因 特异 性 dsRNA 后 发 现 , 灰 飞 乱 翅 出 现 残 
缺 ,前 后 翅 不 能 正常 伸展 , 翅 面 不 能 正常 黑 化 。 在 本 
实验 中 EcR 在 有 翅 成 蚜 中 表达 量 也 显著 高 于 无 翅 
成 是 和 若是 ,所 以 我 们 推测 Burs 和 EcR 对 调控 烟 蚜 
翅 展 及 表皮 着 色 有 一 定 的 关联 ;EcR 对 烟 是 翅 的 伸 
展 发 育 是 否 具有 直接 调控 作用 ,是 否 与 Burs 有 关 ? 
我 们 后 续 将 进行 相关 实验 的 验证 。 

本 研究 利用 高 通 量 测序 技术 对 烟 蚜 若 蚜 .无 翅 
成 蚜 和 有 翅 成 是 的 转录 组 进行 测序 和 分 析 , 揭示 了 
烟 蚜 发 育成 有 翅 成 蚜 时 体内 有 关 翅 发 育 的 基因 的 整 
体 表 达 特 征 。 初 步 研究 发 现 , 在 烟 蚜 翅 的 发 育 过 程 
中 ,关于 好 的 生长 与 发 育 的 基因 (又 化 激素 基因 ) 表 
达 发 生 了 明显 的 变化 。 该 结果 为 深入 人 研究 烟 蚜 怒 型 
分 化 机 制 以 及 相关 基因 功能 及 调控 等 英 定 了 一 定 的 
基础 。 
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